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-^ ' In the inflationary cosmology it occurs frequently that the inflaton field is trapped 

in a local, transient minimum with non-zero vacuum energy. The difficulty regarding 
^ \ the curvature perturbation produced during such a stage is that classically the inflaton 

K*" \ does not move so that the comoving hypersurfaces are not well defined at linear order 

in the scalar field perturbation. In this paper, assuming a mechanism of trapping which 
^+ , resembles a high temperature correction to the potential, we explicitly calculate for the 

■^ \ first time the resulting power spectrum of the curvature perturbation by evaluating the 

t:J" ' quantum two-point correlation function directly. The spectrum is steeply blue with the 

^ . spectral index n-ji = 4. 
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1 Introduction 

Now it is widely believed that inflation [1, 2] takes place at the earliest moments in the history 
of the universe and that after inflation the initial conditions are all satisfied necessary for the 
successful hot big bang universe. One of the greatest triumphs of inflation is that we can natu- 
rally derive a nearly scale invariant spectrum V-ji of the comoving curvature perturbation TZc [3] 
which is required by the recent observations, including the Wilkinson Microwave Anisotropy 
Probe 5-year data [4] where n-n ~ 0.96. This nearly scale invariant spectrum is generated under 
the slow-roll approximation where the infiaton field is very slowly evolving towards the global 
minimum of its effective potential. The calculation of V-n is now a well established subject [5, 6] 
and has become refined with high accuracy [7]. 

However, the slow-roll phase is not a necessary condition for inflation and thus the infla- 
tionary prediction of a nearly scale invariant V-ji is not necessarily true. A typical situation 
where is not slowly rolling off the potential is when it is trapped in a local minimum, i.e. false 
vacuum: for example in the original scenario of inflation [1] it is assumed that is confined in 
a local minimum and the inflationary epoch ends by quantum tunneling. Moreover, although 
this original scenario is observationally not viable, such a phase needs not be ruled out a priori. 
It may lie far outside the observable regime of the last 60 e-folds of expansion out of the whole 
period of inflation-'-. Also, in the case of thermal inflation [9] a short period of inflation is 
provided by a constant vacuum energy due to a temperature effect at the end of conventional 
inflation. 

A problematic fact is that, for such a period the standard calculation of TZc does not work. 
The reason is that when is trapped in a false vacuum, classically = 0^, so that the comoving 
curvature perturbation, which is given by 

7^c~f<50, (1) 

where H = a/a is the Hubble parameter, is not defined. The form of V-r, has been roughly 
guessed [10] but its exact functional form has not yet been known^. But this never means that 
the situation itself is singular, but that we need to adopt a different way of calculation to obtain 
Vtz : we should derive the final result without resorting to the classical homogeneous scalar field 
background [12]. 

In this paper, for the first time we explicitly calculate Vn and the corresponding spectral 
index, n-fi, from a stage of false vacuum inflation. What is important is that the perfect de 
Sitter phase does not last forever. It should eventually end. There are a number of ways to 
terminate this pure de Sitter expansion. Here we adopt a mechanism like thermal inflation. To 
be speciflc, we consider an effective mass-squared which consists of a negative constant term 
corresponding to a bare mass-squared term and a positive term proportional to a~^. The latter 



■l-Note that in the so-called locked inflation [8], most of the observationally relevant part of the universe exits 
the horizon when (/) is effectively trapped in a transient local minimum. However, the stable false vacuum is 
supported by the rapid oscillation of another scalar field coupled to ip. This is different from what we are going 
to discuss. 

sAs there is no background evolution of cj), it may not be quite proper to call it 'infiaton'. Nevertheless, we 
call it the infiaton field simply because its potential energy at the false vacuum is the cause of the inflationary 
de Sitter expansion. 

^Note that in Ref. [11] the spectrum of the field fluctuations, P^^, is calculated. 



term is equivalent to a temperature effect g'^T'^ oc a^^ [13], i.e. such a term can arise due to 
possible couplings to thermal bath. 

This paper is outlined as follows. In Section 2, we write the two-point correlation functions 
of the inflaton and the energy density. In Section 3 we first calculate the power spectrum of 
the gauge invariant intrinsic spatial curvature perturbation $ using the two-point correlation 
function calculated in the previous section. Then we extract the final form of Vn- Finally, we 
conclude in Section 4. 

2 Two-point correlation functions 

2.1 Inflaton fleld two-point correlation function 

We consider a theory with the Einstein-scalar Lagrangian, 

L = ^R - \g^-d,(l)d,(l) - \/(0; t) , (2) 



2 2^ 



where mp, = SttG and the potential is assumed to have the form. 



tpi 
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V = Vo+^-ml^4>\ (3) 



where 

with Vq > 0, m?. < and /x^ > 0. We consider the stage when the effective mass-squared is 
positive, m^g > 0, so that the inflaton is classically trapped at = 0. The background Hubble 
parameter is given by 

3iJ^ = ^ , (5) 

and the cosmic scale factor during this stage can be well approximated by the pure de Sitter 
expression, 

a = a^,exp[iJ(t - t^)] = — -— , (6) 

— Hti 

where t* is an arbitrary fiducial time and 7] is the conformal time. Note that a condition for 
= to be sufficiently stable is ml^/ H"^ ^ 1, hence we must have 
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^,y » 1 . (7) 



In the following we focus on this stage. We also note that it is generally assumed that \mV\ > H^ 
in the case of thermal inflation. 



Since is trapped in a false vacuum so that (0) = during this stage of pure de Sitter 
expansion, the scalar field perturbation is in fact equal to the scalar field itself, i.e. 

50 = 0- (0) = 0. (8) 

Now we begin with considering the two-point function 

G{x,x) = {(f){x)(f){x')) 

= i(0(x)0(a;') + 0(x')0(a:)) 

= Ig^'\x,x'), (9) 

where G^^\x, x') is the symmetric two-point function. We have an exact expression for G'^^^x, x') 
as [14] 

rWl ', ^' r, u ^ l+P(2coshg-2Z)V^ 

G^'(x,x)= — 7T / ascoshius) — -. ; ^,„ ,„ — 

^ ' ^ 27r2 Jo ^ ^ (2coshs-2Z)3/2 

X exp [-p{2 cosh s - 2ZY^^] , (10) 

P=\fWm', (11) 

Z=t±£^jl^ (12) 

r'^=\x-x'\^, (13) 

Note that the form of G*^^-'(a;, x') of our interest is the one in the limit p ^ 1, or equivalently in 
the early stage of inflation r] -^ — oo", when is trapped in a transient local minimum. 

A technically important point is the existence of the term exp [— p(2coshs — 2Z)^/^] in 
Eq. (10). Since cosh a; is exponentially increasing as x increases and we are interested in the 
limit p ^ 1, this term is highly suppressed for large s, making contribution to the integral from 
this region negligible. Thus, the dominant contribution of the integral comes from the region 
near s ~ 0, and hence we can expand the hyperbolic cosine function around this region and 
take only the leading term. Using coshx = 1 + x'^/2 + ■ ■ ■ , we have 



where 



and 



G'^^\x,x') 



l+p[2(l + sV2)-2Z] 




1/2 



[2(1 + sV2) - 2Z]3/2 

1/2- 

(15) 



"Note that the opposite hmit 77 ^- is discussed in Ref. [12] and the spectrum is consistent with the standard 
resuh Vn = [HV(27r0)]2. 



Another point to be kept in mind is that we are uhimately interested only in super-horizon 
scales, just as in the case of standard slow- roll inflation. That is, the two points x and x' are 
space-like separated with their distance being much larger than H^^. Thus introducing a new 
variable 



u^l-Z = 1' , " . (16) 



we find for m ^ 1 



2 



TTZ poo 

G('\x,x') ^ f-f / dsexp [-Pis' + 2uY/'] 
H' p 



2 K,{pV2u), (17) 

27r^ ^/2u 

where we have used an identity of the modified Bessel function of the second kind Ki,{x), 



Now, using the asymptotic form 



^-xVWT^^^ = zKi{xz) . 





2;>i V 2z 



we have 



^ ' ' 27r2 V 2 (2m)3/4 ^ ' 

Therefore, the two-point function in the regime of our interest is given by 

G{x,x') = -G'^^\x,x') 

2.2 Energy density two-point correlation function 

Although there is no classically evolving background scalar field at the stage of our interest, 
this does not mean that there is no energy density fluctuations. In fact, because of the quantum 
vacuum fluctuations of the infiaton field, there exist fluctuations in its energy- momentum tensor. 
To evaluate the curvature perturbation from this stage, calculating the energy-momentum 
tensor in the pure de Sitter background is not sufficient. We have to take into account the 
metric perturbation. But under the situation of our interest where there is no background 
evolution of 0, it is exactly the same as the standard quantum field theory in curved space- 
time. Thus there is no metric perturbation at linear order in the field fiuctuations 50, or the 
scalar field itself 0: see Eq. (8). So in the previous section it is perfectly legitimate to consider 
in the given homogeneous and isotropic background. The metric perturbation 6g^i, appears 
at second order in 0, i.e. it is linear in the perturbation of the energy-momentum tensor. Then 
we can apply the standard linear perturbation theory. This is what we are going to do in this 
section. 



The energy-momentum tensor is given by 

T,, = d,(t>d,(t> - Ig^u (^"^9,09^0 + 2V) , (22) 

and we set 

(5T^ = T^ - (T^) . (23) 

It is important to note that, as can be read from Eqs. (22) and (23), the hnearly perturbed 
energy-momentum tensor ST^iy is quadratic in 0, at which order the metric perturbation comes 
into play and we can follow the standard cosmological perturbation theory: the gauge-invariant 
density perturbation A, which is the density perturbation on the comoving hypersurface on 
which T^i = 0, is expressed as [5, 12] 

V\pA) = V\-T%) + ?,Hd\-T\) , (24) 

where V^ = 6^^didj and we have chosen the time coordinate to be the cosmic proper time; 
x° = t. 

Now we introduce the two-point correlation function of V^(pA), 

D{x,x') ^ {Vl[pA{x)]Vl,[pA{x')]) 

= frit)fr'''^'md^' {[d,d^'dp>G{x,x')] [d^>d,dM^,x')] 



+ [d,d,,G{x,x')] [d^d,d^,dp,Gix,x')]} . (25) 



fOjO _^xj 

Ji 2 i 


J 


ir = - 5/ , 


ft' =a-' 


S^S'^i + i (5>5^-' + 6,^6^'^) 



Using the expression for the energy- momentum tensor, Eq. (22), the coefficients ff'^" are found 
to be [12] 

/r' = f-'' =i;S/ , (26) 

(27) 

(28) 

and zero otherwise, where the potential dependence is eliminated in favor of the spacetime 
derivatives using the field equation for the scalar field. Then substituting the coefficients into 
Eq. (25), collecting non-zero components, expanding the Kronecker delta terms and finally 
rearranging the indices, we obtain a rather lengthy expression given in Appendix B, Eq. (84). 
There the time derivatives are those with respect to t. For later purpose, it is convenient to 
express the time dependence in terms of the conformal time r] defined by dr] = dt/a. Since the 
two-point function depends only on the comoving distance between the two points r = \x — x'\, 
we may then express the coordinate dependence of D as 

D = D{r;r],r]') = D{r;r]',r]), (29) 

where the second equality comes from the fact that D is symmetric under the exchange of r] 
and f]'. 
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Now, inserting Eq. (21) into Eq. (84), we find after some amount of calculations that, inter- 
estingly, the most significant contribution comes from the terms without any time derivatives, 
i.e. the terms multiplied by a~*^ in Eq. (84). To leading order, the two-point correlation function 
D(a;, x') evaluated at an equal time r] = rj' is given by 

/^(a:,a:')-(^)'l6Mi/#|^A-^-. (30) 

We note that we can obtain the same result by substituting the exact expression Eq. (10) into 
Eq. (84) first and then making use of Eqs. (87) and (88), and finally collecting the leading 
terms. 

3 Power spectra 

3.1 Power spectrum of $fc 

With the two-point correlation function given as Eq. (30), we next turn to its Fourier transfor- 
mation. Denoting the Fourier transformation of a function /(r) to be 

mik)^lj^J{r)e~^'-, (31) 



then 



^^^.-,„„,^(m)^8 /■ <'''- «""■'-./..■. 



— 167r(ifr7)4^/-— ^ / dr jo(A:r) , (32) 



yiTi) ' " /i3 ^ (27r)3/2 y^ r 

where the spherical Bessel function jo is given by 



jo(a:) = . (33) 



X 



Notice that the function exp{—2fir)JQ{kr)/r blows up to infinity at r = so the integral does 
not converge. 

An important point to remember at this stage is that, we are interested in the correlations 
of two points which are separated by super-horizon scales, r ^ \ri\. Thus, the singularity at 
r = should not matter and we may introduce a cutoff at a small r for the range of integration. 
Since n'^r]'^ ^ 1 by assumption. This implies that the region of our interest satisfies /xr ^ 1. 
Hence, a natural choice of this cutoff scale would be l//i. Also, since we are interested in very 
large scales, i.e. very small k regions, we can expand sin(A;r) = kr — {kr)^/3\ + ■ ■ ■ . Then, with 



the modified integration range, we have 



dr Jo{kr) ^ / dr Jo{kr) =— / dr — -— sin(A;r) 

^ Ji/f, r kJi/^ H 

1 r , e-^M'- / Pr3 

dr — —- kr 



k Ji/^ r2 V 6 



1 (k 



2 



Ei(-2)- — (-) , (34) 



where we have used the definition of the exponential integral function 

-Ei(-x) = / rft = -7 - log X - V ^^ , (35) 

-'^ n=l 

with 7 PS 0.577216 being the Euler-Mascheroni constant, so -Ei(-2) k, 0.0489005 > 0. It 
should be noted that the argument —2 here is due to our choice of the lower cutoff of the 
integration, r = l//i. In general one may choose any value for the lower cutoff of r as long 
as r = (9(l//i). Hence we should not regard the actual value of — Ei(— 2) to be quantitatively 
meaningful. Instead we should regard it as giving a factor of order unity. In any case, as clear 
from the above, if we are interested in the range of k such that k <^ n, that is, the modes 
which leave the horizon before the time —r] = 1/fi, the second term proportional to (k/fi)'^ is 
negligible. 

Now we can explicitly write the power spectrum of $, the gauge invariant intrinsic spatial 
curvature perturbation in the Newtonian (or longitudinal) gauge. From the perturbed Einstein 
equations, we have the well known relation [5] 

$ = -4nGpA , (36) 

where the factor pA shows up which appears in the definition of D[x,x'), Eq. (25). Hence 

pA = -T^^$ = -2ml,{Hr^rv'<!> . (37) 

Since we are interested in the correlation function of two points which are apart on super- 
horizon scales, the leading contribution of the spatial gradient on the function of the form 
Q-'if^r i^z gives 

'g-2Mr\ I Q 



V 



Thus, substituting Eq. (37) into Eq. (25), 



dr \ r-^ 



g-2^r 

{2py—^ . (3j 






Dix, x') = Aml,{Hr^Y (V^ [Vl^{x)\ Vl [Vl<^{x')\ ) 

^ 4m^i(ffr7)^(2/i)«($(x)$(x')) . (39) 



Therefore, equating Eq. (30) with Eq. (39), the two-point correlation function in configuration 
space is given by 

^^{r) =mx)<^{x + r)) 



4m|.i(iJr7)4(2^)8 \2n J ' " (/xr) 

2tir 



^ [ H vm' 



64 \27rmpi/ (/ir)^ 

where we can see that it is exponentially suppressed. Now taking the inverse Fourier transfor- 
mation, we can relate ^$ and P$ as 

^* = |^y"rfVe*(r)e-^'»•^ (41) 

which includes the integral J„ e~'^'^^JQ{kr)/r that we have already calculated and is given by 
Eq. (34). We have thus 

/.3.,,,,-....._.,_,^(_JL_)*M!. „, 

Note that we have fi\ri\ ^ 1 so that the temperature term, the second term of Eq. (4), dominates 
the effective mass-squared, and hence $ oc r^^. This implies that $ we are calculating is decaying 
in time, hence it may seem that it does not contribute to the final power spectrum at all. 
However, it turns out that this behavior of $ correctly corresponds to the growing adiabatic 
mode as we shall see later. To summarize, we have the power spectrum of $ to leading order 
as 

To translate this into the power spectrum of the comoving curvature perturbation V-n, we need 
to know (p + p) which is classically 0. 

3.2 Power spectrum of the comoving curvature perturbation 

3.2.1 (p + p) 

As it is well known, the adiabatic density perturbations responsible for the large scale structure 
of the universe today is represented by the curvature perturbation on comoving hypersurfaces 
TZc on super-horizon scales. However, the definition of the comoving hypersurface, T^i = 0, 
becomes meaningless in the pure de Sitter space. This is simply because T^i, oc gfj,u in pure 
de Sitter space, hence T'^i is identically zero. In other words, there exists no preferred rest 
frame in pure de Sitter space. Nevertheless, in the present case, we do have a preferred frame 
because of the time dependence of mg^. Therefore, the vacuum expectation value of T^i, will 
no longer be de Sitter invariant. In particular, we expect p + p to have a small but non- 
vanishing vacuum expectation value. This fact enables us define the comoving hypersurface 
and hence the spectrum of the comoving curvature perturbation, Vn- Here we assume that the 



thermal contribution to {p + P), which is proportional to T'^, is small compared to the vacuum 
contribution from 0. We will come back to this point at the end of this subsection. 
From the energy-momentum tensor for a scalar field, we find 



.0 li2 , l(V0f 
2^ 2 a2 



P = n= 00+0^^ + ^0), (44) 

3*2^ 6 a2 



P=-in = ;0^-ii^-n0). (45) 



Thus we have 

iW)2 

3 a 
Using the standard Fourier mode expansion of the scalar field, 



P + P = ^' + 7;'-^- (46) 



X 



"^'^ afc0fc(77)e*^ + al,0*(77)e-*^ , (47) 



(27r)3/2 

we then obtain 

(^=>-(t?>-^/(i5«'M/W. (48) 

The explicit form of the mode function (pi^ and the calculation of {(fP) are given in Appendix A. 
We are left with 



')-4^^(!-^'0(ipr''"^^"'«^"'<^ 



vr , , /I 3 ^ 47r '""^ 

\3 



+ 5?<-"' (2 + S;j (^X dkkHe + ,^)HiU^)H':i,u) 



+ £1'-"' (i - I;) (l)J /°° **''*' + A^)C>,(.)hSi(^) . (49) 

In the same way, we find 

1 (V|)^ = ^(_^) J!!_ rdkk'Hi'\z)Hi'\z) . (50) 

Now we have all the ingredients to calculate {p + p)- 

As clear from the above expressions, the momentum integrals diverge at large k. This is the 
standard ultraviolet divergence that should be regularized or renormalized. To regularze the 
divergence, we introduce a simple cutoff at a large physical momentum, k/a < {k/a)c = HA 
where A is a large number. 

Now let us evaluate the integrals in the limit rj -^ —00, i.e. during the early stage of 
inflation. This is the region of our interest. In this case, we have 



^/WT]^{-r]) > p{-r]) :^ 1 , (51) 

so the Hankel functions are approximated as 

/j(i-2)(^) ^ /Ae±<-f-f ) (^ » 1) , (52) 

V TfZ 



thus 

Hi'\z)Hi'\z) r^ - (^»1), (53) 

nz 

which is independent of the parameter z/. Then after some calculations, we find 



H^ A4 /9 2 



I^^f^\rfr_fj^\r]f\^fj^\rf^^ f fi\v\^ ^ 



2 4 V 2A2 y 2 V 2A ^ 

' ' (54) 

To renormalize the above expression, we assume that the de Sitter invariance will be (formally) 
unbroken if /x = 0, that is, if there is no preferred frame. Then the renormalized expectation 
value is given by 

(p + p)rc„ = lim [{p + p) (A, fi)-{p + p) (A, 0)] , (55) 

A>1 

where (p + p){A,p) is the vacuum expectation value given by Eq. (54). This leads to 



{p + P)ren = Y^J^ipV) i 1 + lo? 



2 A/ 



(56) 



This should be valid for 1 ^ fi\ri\ <^ A. 

Now, let us check if the above result is indeed greater than the thermal contribution. As- 
suming there are N effective massless degrees of freedom that are thermal and that couple to 
the (h field, we have an estimate 



■■5 

,2 



(p + P)t ~ NT^ and ^ = g'NT' , (57) 

where the suffix T stands for thermal, and we have assumed the same coupling constant g^ 
for all the thermal fields for simplicity. On the other hand, the contribution from the vacuum 
fiuctuations of 0, Eq. (56), may be approximately expressed as 



p 



2 



(p + p)^~|mj| ^±^=\ml\g^NT\ (5J 



Hence in order for the contribution to dominate, we must have T^ <^ |m^| g"^. Since the stage 
of our interest is |m|| < fx'^/a'^ = g'^NT'^, this leads to the condition, 

T^^\ml\ g^ <g^NT\ (59) 

This can be satisfied only if we have g^N ^ 1. Namely, the dominance of the contribution 
over thermal corrections can be realized if there are sufficiently large massless degrees of freedom 
which are in thermal equilibrium at an early stage of the false vacuum infiation. However, we 
should mention that the condition g'^N ^ 1 implies that the theory is in a regime of strong 
coupling. Thus the thermal contribution can dominate in general, suggesting that thermal 
corrections may be worth investigating in more detail, although it is outside of the scope of the 
present paper. In the following subsection we assume that thermal corrections are negligible. 
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3.2.2 Extraction of growing solution 

Using the results obtained until now, we can compute the comoving curvature perturbation TZc- 
First let us recapitulate the standard result. 

On super-horizon scales, for a completely general equation of state, we have the general 
solution for $, 

$ = ttC*!— / (1 + w)a^{7]')d7]' + decaying mode , (60) 

where Ti = a' /a, w = p/p and the decaying mode is proportional to H/a^. The initial time rji 
can be chosen arbitrarily since its change only affects the contribution of the decaying mode. 
Given the Newtonian curvature perturbation $, TZ^ is expressed in terms of $ as 

2t' + (5 + 3.)«t (61) 

Plugging the general solution given by Eq. (60) into the above, we find, neglecting the decaying 
mode part, 

7^e = Ci . (62) 

Thus the coefficient Ci indeed corresponds to the amplitude of the growing adiabatic mode. 

In order to extract out the final curvature perturbation amplitude Ci from our calculation, 
we need to perform the integral in Eq. (60) explicitly. To do so, we need the information of a 
and 1 + w. At leading order approximation, we may assume the universe is de Sitter. Hence 
we may approximate the scale factor and the Hubble parameter as 

^ n = —. (63) 



—Ht] —1] 

As for 1 + w, however, we must take into account the small deviations from de Sitter. From 
Eq. (56), we note that the logarithm log[(/i?7/2A)^] is negative for fj,\ri\ <^ A. So, neglecting the 
slow logarithmic behavior and take into account the fact m^ < 0, we effectively have 

(p + p)ren = AY^^(/ir7)^ (64) 



where A is 0{1) and positive. This gives 



A \mV 



^+'^' = wi^<'"''- '«' 



This should be valid for fj, r] ^ 1 as long as 1 + w < 0(1), that is. 



m 



2 

PI ^-^ ..2 2 



I 2 I 

\mi\ 



> pY > 1 • (66) 



Keeping in mind this range, the integral on the right-hand side of Eq. (60) can be easily 
evaluated to give 

2 A l"^il 

$ = —CiK^pr]) + decaying mode ; k = ^ — — . (67) 
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Tn(k) = ^:S^. (6i 



Now it is easy to obtain the final amplitude of TZc- We just have to divide the early time 
solution for $ by K,{firi)'^. Then the spectrum of the conserved comoving curvature perturbation 
is 

Inserting Eq. (43) into the above, we thus finally find the spectrum of the comoving curvature 
perturbation at A; ^ ;U as 

where B = — 2Ei(— 2) ^ 0.0978010. Thus, with Vnik) oc A;"^^^, where n-pi is the spectral index, 
we have 

n7e = 4, (70) 

i.e. a very blue spectrum. This result is in accordance with a naive expectation. That is, 
for a very small k which leaves the horizon when = is strongly stable, the generation of 
the curvature perturbation must be severely suppressed. On the other hand, as k approaches 
/i, the amplitude of the curvature perturbation should increase. This is because the effect of 
the negative m| gradually becomes more and more important as the universe expands and the 
instability sets in at 

I 2 I 
2 2^ 22 \"''4>\ /«-, \ 

^T] <fir]^ = -j^ . (71) 



c- 



It is then expected that the fiuctuations would become very large dX r] = 7], 

4 Conclusions 

In this paper, we have calculated the power spectrum V-ji and the corresponding spectral index 
n-ji of the comoving curvature perturbation TZc produced when the infiaton field is trapped 
in a local minimum of the effective potential with non-zero vacuum energy. The difficulty we 
confront when we adopt the conventional approach of the calculation of TZc is that classically the 
infiaton is well anchored if |?72eflf | ^ 3if/2 so that = 0, and accordingly comoving hypersurfaces 
on which TZc is given become singular: this is because we cannot define comoving hypersurfaces 
when = 0. 

To evade this difficulty, we have used a pure quantum field theory approach to calculate the 
two-point correlation function of the infiaton field and the perturbation of the energy density 
in the exact de Sitter background. This de Sitter phase, due to a non-vanishing false vacuum 
energy, shoud not last forever but should eventually be terminated so that the standard hot 
big bang evolution of the universe can commence. We achieve this by adding a comoving mass 
term, or equivalently a thermal correction, ^ jd^ ^ to the potential as shown in Eq. (4). This 
breaks the perfect de Sitter invariance, and allows the vacuum expectation value of p + p to 
be non-vanishing. Then given the fact that (p + p) 7^ 0, we have explicitly calculated the final 
comoving curvature power spectrum Vn^ given by Eq. (69). The spectral index is found to be 
very blue, n-ji = 4. 
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We believe our results are widely applicable: for example, we can directly constrain the 
production of the primordial black holes after thermal inflation [15], which has not yet been 
studied anywhere including the original references [9]. Also since the curvature perturbation 
is quadratic in the scalar field it is highly non-Gaussian. This issue will also be reported 
separately. 
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A Inflaton field mode functions 

Since there is no classically evolving background, the quantization of the inflaton can be done 
without worrying about the metric perturbation, that is, as in the standard quantization of a 
scalar field in curved spacetime. 

We start with the Fourier expansion of the infiaton field. 



[X 



au<PMe'''-^ + alcl^^Me-^'^-^ , (72) 



(27r)3/2 L 

where the annihilation and creation operators ap, and a\, satisfy the standard commutation 
relation 

[afc,at] =5'^^\k-q). (73) 

The mode function 0^ is called the positive frequency function. It determines the vacuum 
annihilated by a^- The mode equation is 



Ul 






{e + ^?) + 



2 -\ 

ml 



jy2^2 



0, (74) 



with the normalization 



0fc0fc - 4>*k4>k = -\ = -tH'v' , (75) 



where a prime denotes a derivative with respect to rj. 

The positive frequency function 0^ appropriate for the infiationary universe is given by 



MV) = ^Hi-r^f'Hl'^ (_^/pT7^^) , (76) 
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where v is given by Eq. (14), 



'9_m^ / 3 
iJ2 Y 2 



(77) 



and Hi, (2;) is the Hankel function of the first kind. This mode function satisfies the asymptotic 
boundary condition that it reduces to the one for Minkowski vacuum in the hmit ri -^ —00, 
corresponding to the high frequency hmit where the cosmic expansion can be totaUy neglected. 
Below, we list a few formulas for 0^ and its derivatives which are used for the computation of 
the components of the energy-momentum tensor. 
Taking a time derivative of Eq. (76), we find 



^'^(^) - T, 



^H{-^fl^H^^ {-y^WTJT'V 



Hi-V) 



1/2 



-^Hi'^ (^-,/WT7v) - Vk^TJi^i-v)Hi'^' {-V^T7 



(75 



where the prime of the Hankel function means the derivative with respect to the whole argu- 
ment, i.e. Hu (z) = dHp \z)/dz. 

Using Hiy (z) = Hu {z) for real v and z, we find 






where z = ykP^+fj? [—r]). Eliminating Hl{z) by using the Hankel function identities 

H,_^{z) - H,+,{z) 



H'A^) 



V 



r(i) 



which hold for both Hi, iz) and H{, (2;), we obtain 



H,{z)=^[H,^,{z) + H,+,{z)] 

(2), 



(79) 

(80) 
(81) 
(82) 



n 



MvHl'iv) =-^H\-r^) 






^3) 



B Formulas for energy density two-point function 

Here we give an explicit expression for the two-point function D{x,x') introduced in Eq. (25) 
in terms of the scalar field two-point function G{x,x'). We also give useful formulas for the 
spatial and time derivatives when they act on G{x,x'). 
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Substituting the expressions for the coefficients fj^'^ given by Eqs. (26), (27) and (28) into 
Eq. (25), we find 

D{x,x') = d'd'' mdodo^G) {d.^dodo'G) + {dodo'G) {d^d^^dodo'C)] 

- (So^G) (d.dydoG) - (dydoG) {d,dod',,G) 

- {did,,do'G) (do'd^oG) - id.do'G) {dydo'd^oG) 
+ (a,4G) (dyd^oG) + {^^^,,G) (d^odo'G)] 



+ a 



-2 



d'd^ 



d' dk'doG (d^dydoG) + id,>doG) d,d' dk'doG 



+ (dyd'^'doG) (dA'doG) + (dk'doG) (d^dyd'^'doG 



+ a 



d'& 



d^d'^'d^.G) {dydlG) - [d,d.yG) [d'^'dk.dlG 
d.dyd'^'G^ [d^^dlG] - [d,dk<G) [dyd'^'dlG 

+ [d,d,.d^,G) [d^d^'d^.G) + (S.So'G) {d^d''d,.d^,G) 
+ {d^dyd^,G) {d,d''d^,G) + {dkd^.G) {d,d''d,,d^,G) 

- {d,dl,G) [d^d^dyG) - [d^dyG) [d^d^dlG) 

- [d^dlG) [d^d'^d^.G] - {d^dyG) (9,9^92 G)] } 
(d,dvd''G^ [dkd'^d.^G) + {d,d,,G) (d^d'^dyd^'G) 

+ {didyd^'G\ [dkd^dvG] + {d^dyG) (d^d^d^.d^' G 

+ [d^dyd^'G^ (9,9"9,vG) + [d^dyG) [d^d^'dyd^'G 

+ [dkdyd^'G\ [d.d^dvG) + [d^dyG) [d.d^dyd^'G 



^4) 



Note that we have not used any special properties of G(x, x') in the above. 

To proceed, we first make use of the properties of G(x, x'). As far as the spatial dependence 
is concerned, it depends only on the distance between two points r = |a; — a3'|. Hence we have 



•9^' 



-d^. . 



^5) 



As for the time dependence, G(x, x') is symmetric in the interchange of t and t' . Thus we have 

G = G[r;t,i!) = G[r;i!,t). (86) 

Now, for the two-point function G(x^x'^ with the above properties, the spatial derivatives 
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acting on G{x, x') may be expressed as 



di = —dr , 
r 

rf '' ryj rf "^ / rj-i '^ rf '' /yj -i /y /yj /y \ / I I 

didjdk = :^—&r + ^ij ^ ^jk— + Ski ^ ] {-^r ^^r I , 

f '' lyj ly "^ ly •' 



didjdkdi = "^ -"j^ •" g4 



ly" lyJ ly'" /y' 

lij lij lij JU 

r 

Xtxi ,-yit lyJ ly'' /yj /y'^ ly '' /y" /y'' /y'^ /y'' lyj t-v /y'' lyj ry'^ /y'' \ T 

+ I dij — — + Oik — 5- + Oil — 5- + Ojk^ir- + C)j7 — r- + Oki — 5 T -a^ 



+ 



rf ly lyJ ly lyJ ly ly ly ly ly 



SijSki + Sjk5ii + 5ki5ji - 3 5ij — — + 5ik — ^ + Su — ^ + ^jk^ir + ^ji — 5- + ^k 






'It r, I "Jfc r, I "J( r, I "ifc r, I "l( r, I 'J kl r, 

-' rf^ y ^ yZi yZi y^ y ^ 

^A - ^^^^ . (87) 



Similarly, for the time derivatives acting on G{x,x'), using drj = dt/a, we have 
do = H{-ri)dr, , 

9o(9o' = H'^{r]r]')dndr,' , 

dldl, = H\r^i) [(W)5J5J - {-n)dld,, - {-n')d,dl, + d,d,] . (88) 

These formulas are used in the explicit evaluation of D{x,x'). 
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